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1. Introduction 



Supersymmetric gauge theories in four dimensions possess a very rich physical and 
mathematical structure. A surprising recent example is the connection between twisted 
M = 4 theories and the Langlands program [Q. In the 90' s, the Seiberg-Witten solution 
of M = 2 theories and their realization of confinement after breaking to M = 1 0|| was a 
major breakthrough in showing the control supersymmetry gives over quantum corrections. 
Seiberg dualities [|J in M = 1 theories show the power of holomorphy in theories of more 
phenomenological interest. Also, deep connections, such as the AdS/CFT correspondence 
0, have been found between gauge theories and string theories. 

In this paper we explore an unexpected connection between Seiberg-Witten theory 
and the theory of quadratic differentials on Riemann surfaces. The quadratic differentials 
that will be of interest to us are called Strebel differentials 0. Given a Riemann surface 
with n marked points, a Strebel differential induces on the Riemann surface a metric that 
makes it look like n semi-infinite cylinders glued along a graph. This metric is unique up 
to a choice of n real numbers that correspond to the radii of the cylinders. 

Strebel differentials naturally establish a bijection between M. g , n x IR+ 5 the moduli 
space of Riemann surfaces of genus g with n punctures, with a positive real number as- 
sociated to each puncture, and the moduli space of metric ribbon graphs - ribbon graphs 
with a length associated to each edge [|7|||. Because of this property, Strebel differentials 
have been useful in string field theory [j^PIPOfl and have played a central role in the cal- 



culation of tachyon amplitudes [ |IT| , |T!^ , |T^ , |T4l , |To| . More recently, they have appeared in a 
proposal by Gopakumar [pT6| , p77| , p~8|] to relate free field theory ribbon diagrams to closed 
string worldsheet correlators. See also [ fIT)| , |2Ti| , pT| . 

A common roadblock in the different attempts to use Strebel differentials in physics is 
the difficulty in constructing such differentials explicitly. This is because the condition for 
a quadratic differential to be Strebel is transcendental and expressed in terms of elliptic 
functions. Very few explicit differentials are known, although some numerical []T3|JT5f| and 
perturbative approaches |2D| have been developed. 

In this paper we restrict our study to Strebel differentials with some integral prop- 
erties. More explicitly, we impose that the lengths of the edges of the associated ribbon 
graphs be integers. This condition leads to a huge simplification of the problem: it replaces 
transcendental equations by polynomial equations. The idea of restricting to Strebel dif- 
ferentials with this property came from physics; in particular, from supersymmetric gauge 

1 



theories, where a generating function of chiral operators can naturally be interpreted as 
an abelian differential on a Riemann surface with integral periods f22] , p3"l , |24|| . 

At first it might appear that imposing the integrality of the lengths is a very strong 
constraint and that such differentials are very scarce. This, however, is not the case and 
in later sections we show that one can approximate any Strebel differential, to any desired 
degree of accuracy, by one that can be obtained by solving polynomial equations. 

Moreover, we will show that constructing Strebel differentials with integer lengths is 
equivalent to solving for a particular class of Seiberg-Witten curves with Argyres-Douglas 
singularities [^5| in the moduli space of a U(N) M = 2 supersymmetric gauge theory with 
massive flavours. More explicitly, we replace the transcendental equations by the following 
factorization problem 

y 2 sw = P 2 (z) + B(z) = Q\z)R\z) (1.1) 

where P(z), B(z), Q(z) and R(z) are all polynomials. Although much simpler than the 
original transcendental equations, solving ( |1 . 1| ) is still a challenging problem in general. 
However, it turns out that the relation to Strebel differentials allows us to write a differ- 
ential equation that reduces the problem, by means of solving linear equations, to only a 
small set of polynomial equations. 

In the case we mostly study, which is the case with four punctures, the number of 
polynomial equations is always two. These are solved by computing the resultant of the two 
polynomials. These resultants have interesting factorization properties over Q, connected 
to the fact that Strebel differentials with integral lengths can also be thought of as the pull 
back of a meromorphic differential on the sphere by a Belyi map |2^,^| . We will comment 
only briefly on this point, leaving a detailed discussion to a forthcoming publication |2^ . 



This paper is organized as follows. In section 2, we review the definition of Strebel 
differentials, their essential properties and the equations that need to be solved in order 
to construct them. As mentioned above, these equations are in general transcendental. 
We use the relation between Strebel differentials and ribbon graphs to explain how to 
relate Strebel differentials with integer lengths to those with rational lengths. From this 
connection we explain how to approximate any Strebel differential with arbitrary real 
lengths by constructing a related differential with integer lengths. In section 3, we show 
how the problem of finding Strebel differentials with integer lengths is purely algebraic 
and we relate this problem to the factorization of Seiberg-Witten curves with Argyres- 
Douglas singularities. In section 4, we show how to use the analytic structure of a Strebel 
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differential to construct a differential equation, which is the main tool that will allow 
us to solve the examples treated in Section 5 and in the appendices. In section 5, we 
consider examples that illustrate several points discussed in the previous sections. The 
examples considered are the sphere with three and four punctures. In the latter set of 
examples, after discussing the general problem, we specialize to that of equal residues 
(n, n, n, n) with n G N. In section 6, we discuss how the relation to Strebel differentials 
makes its appearance directly in gauge theory and in the string theory realizations of it. In 
section 7, we place our analysis in a broader mathematical setting and summarize work in 
progress that generalizes the correspondence discussed here between Seiberg-Witten curves 
and critical graphs of Strebel differentials to Grothendieck's theory of dessin d'enfants or 
"children's drawings" . Some technical details and additional examples are collected in the 
appendices. 

2. Strebel Differentials 

The central figure in this paper is a special class of quadratic differentials on Riemann 
surfaces called Strebel differentials. In order to understand their relation to gauge theories 
we will exploit many of their properties, in particular their uniqueness, some equivalent 
definitions and their relation to metric ribbon graphs. Therefore, in this section, we will 
summarize a few relevant facts in the theory of Strebel differentials. 

Definition Of Strebel differential 

We are only interested in Strebel differentials defined on a Riemann sphere S with 
punctures. A Strebel differential || is a meromorphic quadratic differential with the fol- 
lowing properties: 

/. It is holomorphic on S \{pi, ■ ■ ■ ,p n } , where the p^s are the location of the 
punctures. 

II. It has only double poles located at z = pi for i = 1, . . . , n. 

III. The set of open horizontal curves is of measure zero. 

A horizontal curve is defined as follows. Denote a generic quadratic differential by 4>(z)dz 2 
and denote by 7 a curve on S parameterized by t. Then 7 is a horizontal curve of (p(z)dz 2 
if (/>(7(t))(7'(t)) 2 Gl+ for all t. 
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A theorem by Strebel guarantees that a differential satisfying all three conditions 
exists. Moreover, if the residues at each pole are given and are in IR + then the differential 
is unique. 

The open horizontal trajectories form the critical graph of <p(z)dz 2 : the vertices are 
the zeroes of <p(z)dz 2 and the edges are open horizontal curves between two zeroes. For 
generic locations of the punctures the differential has simple zeroes and the vertices of 
the critical graph are trivalent. Degenerate configurations with vertices of higher valency 
occur at the points in the moduli space where the differential develops higher order zeroes. 

A quadratic differential <p(z) dz 2 naturally induces a metric \<p(z)\dzdz on the Riemann 
surface. This metric is flat almost everywhere around each pole of <p(z)dz 2 , except on 
the zero measure set of open curves. The poles are pushed to infinity and the surface 
looks like a collection of n semi-infinite cylinders glued along the critical graph of the 
Strebel differential, where all the curvature is concentrated. The residues of the differential 
determine the radii of the cylinders. 

A Strebel differential can also be completely characterized by giving all the Strebel 
lengths, i.e. the lengths of the edges of the critical graph measured in the Strebel metric. 
The following condition is equivalent to the third condition above and is the one usually 
implemented in practice to construct a Strebel differential: 

Illb. The lengths a/ 4>{z)dz , with a and b zeroes of the differential, are real. 

Note that this condition makes sense, since we have imposed that the residues of the poles 
be real as well. 

2.1. Finding Strebel Differentials 

In practice, given the location of the punctures {pi,---,p n } an d the residues 
{mi, . . . , m n } one constructs a differential of the form 

Hz)dz 2 = ^dz 2 , (2.1) 

where Q(z) and B(z) are polynomials of degree 2n — 4 and n respectively!!. Usually, one 
uses the £X(2,C) symmetry to fix three of the poles to be located at {0, 1, oo}. Taking 

1 The degree of Q(z) coincides with the number of simple zeroes of a generic Strebel differential. 
In the generic case all vertices are trivalent, so using the appropriate relation between the number 
of edges and vertices in Euler's formula for a planar graph one finds that, if n is the number of 
faces, the number of zeroes must be 2n — 4 . 
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the remaining poles to be at pj, with i = 2, . . . , n — 2, and the roots ( equivalent ly, the 
coefficients) of Q(z) to be unknowns, the Strebel differential takes the form 

47T 2 Z ^ Z -l)2 U]J(z-Pj) 2 

Here, we have chosen the overall coefficient so that the residue at infinity is given by m^ . 
By imposing that the residues, which are calculated as the usual residues of simples poles 



of the abelian differential y/cf)(z)dz, are equal to m; , we can fix n — 1 coefficients of Q{z) 
by solving linear equations. Therefore, n — 3 coefficients are left unknown. Notice that 
y 2 = Q(z) is a Riemann surface of genus n— 3 and therefore it has 2n— 6 independent cycles. 
This is also the number of independent open trajectories of the differential. Imposing that 
these 2n — 6 lengths be real completely fixes the n — 3 complex parameters left in Q(z) 
as functions of the locations of the punctures and the residues. Here is where all the 
complexity of the problem resides. In order to impose that the lengths are real one has to 
find the roots of Q(z), denoted by {zi} and demand that 

Im f f 3 y^(z)dz) = 0. (2.3) 



These equations are in general very hard to solve and this is what motivated us to search 
for simpler equations!!. Interestingly, we will see that imposing the stronger condition of 
the lengths being integral simplifies the problem considerably. 

2.2. Example: Sphere With Three Punctures 

As a first example, let us consider the case with three punctures. This case is simple 
because, since n = 3, after fixing the residues there are no unknown parameters left in 
Q(z). This is clear from the fact that the only length to be computed is given by an 
integral that can be deformed and evaluated as the sum of the residues, which are real by 
definition. Therefore, the lengths are automatically real. If the residues are chosen to be 
{mo, mi, moo}, h can be shown that the strebel differential is completely fixed in terms of 
the residues to be 

md z 2 = -gs p I -h z 2 . (2.4) 

If n > 3 the problem of finding the Strebel differential will also involve imposing conditions 
on the lengths ( |2.3| ) which is much harder since it involves solving transcendental equations. 



A perturbative approach to computing the Strebel differential was initiated in [2C]. 
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2.3. Relation To Metric Ribbon Graphs 

It is important to keep in mind the relation between Strebel differentials and ribbon 
graphs. This relation will allows us to learn some nontrivial properties of the Strebel 
differentials with integer lengths and the polynomial equations that give rise to them. 

The precise statement, specialized to genus zero, is that Strebel theory establishes a 
one-to-one correspondence between Aio tTl x (the decorated moduli space of a sphere 
with n punctures) and the space of metric planar ribbon graphs, i.e. planar ribbon graphs 
with a length associated to each edge. 

First note that the decorated moduli space A4o, n x 1R+ coincides with the moduli space 
of Strebel differentials on the sphere with n punctures, since we can interpret the positive 
real number associated to each puncture as the residue of the differential at that point. 
Strebel's uniqueness theorem guarantees that, for every choice of residues and location of 
the punctures, the differential is unique. Therefore we can rephrase the statement above 
as a one-to-one correspondence between the space of Strebel differentials and the space of 
metric ribbon graphs. 

It is easy to see that the statement is true one way. Given a Strebel differential we 
have defined the critical graph as the set of open horizontal trajectories that connect the 
zeroes of the differential, with the zeroes being the vertices of the graph. This is a genus 
zero graph that has a positive real number, the distance between the zeroes measured by 
the Strebel metric, associated to each edge. This is the metric ribbon graph. 

It is also true that it is possible to associate a unique Strebel differential to any given 
planar ribbon graph, but the proof is more difficult (see 0). 

The critical graph of a Strebel differential contains all information about the Strebel 
differential and the associated Strebel metric. Thinking in terms of the critical graph can 
be helpful to intuitively understand the properties of the differential and we will refer to 
this later in the discussion. In figure 1 we show the critical graph of a Strebel differential 
on a sphere with four punctures that has the topology of a tetrahedron. The associated 
metric ribbon graph can be obtained by thinking about the (thickened) edges as six ribbons 
that meet at four trivalent vertices.! 

3 There are more topologies in this case but we postpone their description to section 6.2. 
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Figure 1: Critical graph of a generic Strebel differential on a sphere with four punctures. 

2.4- Approximating Any Strebel Differential 

As mentioned in the introduction, in this paper we restrict our attention to Strebel 
differentials with integer lengths. Imposing this condition seems to be a strong constraint; 
however, from the relation to ribbon graphs it is clear that this is not so. The condition we 
are imposing on the differential is equivalent to asking that the associated ribbon graph 
has all integer lengths which can be easily satisfied. 

In the rest of this section we give another motivation for constructing Strebel differen- 
tials with integer lengths. We claim that a generic Strebel differential can be approximated 
to any degree of accuracy by a Strebel differential with integer lengths (divided by some 
appropriate integer). 

The basic idea is the following. Consider a generic Strebel differential on a sphere with 
n punctures with an associated ribbon graph with lengths given by {£\, . . . , £3 n _6}.@ Each 
li G M + , therefore we can choose a set of rational numbers that approximate each length to 
a prescribed accuracy. In order words, take e > and let < e , with li e Q. Then let 

us consider the approximating Strebel differential 4>a(z)cIz 2 with lengths . . . ,1^-q}. 
That such a differential exists is clear from the existence of the corresponding ribbon graph. 

4 Recall that the ribbon graph has n faces and that the perimeter of the face is given by 
the residue at the puncture. This is why one usually gives the n residues and {£i, . . . ,£2n-e} 
independent lengths. However, for this section it is best to work with all the 3n — 6 lengths and 
not mention the residues. 
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Let MeNbe the smallest number such that Mli G N for all i. Then <pi(z)dz 2 = 
M 2 <p A (z)dz 2 is a Strebel differential with integer lengths. Therefore, one can construct an 
approximation of the desired differential as 

<j> A {z)dz 2 = -^ ( p I (z)dz 2 . (2.5) 

In the next section we will show how to reduce the problem of finding a Strebel 
differential with integer lengths to a set of algebraic equations and in section 5 we will 
study the series of differentials with residues (n, n, n, n) , with n G N. From the above 
discussion, we can think of this series as giving rise to an infinite set of Strebel differentials 
with residues (1, 1, 1, 1). By explicit computations we will see how the solutions for different 
values of n correspond to a subset of points in a lattice in .M 0,4 • The lattice gets finer as 
we increase n, so that in the limit n — ■> oo, we expect the lattice of points to fill up all of 

-Mo,4- 



3. Strebel Differentials With Integral Lengths 

In this section we show that while finding generic Strebel differentials involves solving 
transcendental equations, if we impose that all the lengths be integers then the problem can 
be mapped to a purely algebraic problem. This condition can also be stated by saying that 



all independent periods and all residues of a/ 4>(z)dz, thought of as an abelian differential 
on y 2 = Q(z), are required to be integers. 

The problem of finding an abelian differential on a hyperelliptic Riemann surface 
with certain kind of singularities and special integrality conditions was encountered in 
H,||,|2|j2|] in the context of U(N) super Yang-Mills theories. More specifically, in ||22| 



the relevant abelian differential was the generating function T{z)dz of (Tr <f> n ) where $ is 
a chiral supermultiplet for a U(N) theory with massive flavors in the presence of a tree 
level superpotential for $ . We follow the same steps as in the construction of |2«| , adapted 
to the particular singularity structure of our problem. In fact, it will turn out that the 
Strebel differential we are after is related very simply to T{z)dz as 



2ni y/(/)(z)dz = -2T(z)dz + R(z)dz . 



for a particular rational function R(z) . 
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3.1. Construction Of The Strebel Differential And Argyres- Douglas Singularities 

We would like an abelian differential on E : y 2 = Q(z) such that all its periods are 
integers. We can achieve this by taking 

y/Wjdz=±d]og(f(z)) (3.1) 
2m 

for some well-defined meromorphic function f(z) on E.i 

We will see below that, given the singularity structure of the Strebel differential (simple 
zeroes and double poles), it is more natural to begin by first constructing a differential of 
the form ( |3.1|) on an "auxiliary" Riemann surface Eq defined by the equation 

y 2 sw = P 2 (z)+B(z) 1 (3.2) 

with P(z) and B(z) polynomials of degrees iV and L respectively. 

From ( |3.1|) it is clear that the zeroes or poles of f(z) of degree k become simple poles 
of ^<p(z) dz with residue k or —A; respectively. Also, since we want (j)(z)dz 2 to be well 
defined on the sphere, for each pole on the upper sheet there must be a pole on the lower 
sheet with the same residue. Since we need poles on both the lower and the upper sheet 
of So , the function we choose to construct the differential (|3.1| ) is 

/(*) = i\ Z \ ~ VSW ■ (3-3) 
1 ! P(z) +y sw K ' 

Every zero of B(z) leads to a simple pole in the abelian differential. Likewise, every zero of 
ysw leads to a zero of the differential. As we will see below, the further requirement that 
the function f(z) be well-defined on E leads to constraints that cut down the number of 
zeroes to exactly the right number for ( |3.1|) to be the square root of a Strebel differential. 

We use an SL(2, C) transformation to set three of the n punctures at {0,1, oo}. 
There are thus n — 1 poles at finite points which we identify with zeroes of B(z), each 
of order m r More explicitly, if {oo, 0, 1,^2, • • • iPn-2} are the locations of the poles and 
{moo, mo, m ii m 2, ■ ■ ■ j are the corresponding positive integer residues, then 

n-2 

B(z) =az mo (z-l) mi Y[(z- Pj ) m i , (3.4) 

J=2 



are e 



Being well-defined on S means that the only monodromies allowed around the cycles of £ 
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where at this point a is an arbitrary constant. With this identification the degree of B(z) 
is L = mo + mi + rri2 + . . . + ra„_2. The degree of P{z), denoted by N, is related to the 
residue at infinity: it is easy to checki that N = (L + m OG )/2 . For convenience, we choose 
the residue at infinity moo to be the largest residue. 

We would like to identify the curve Eq with the Seiberg-Witten curve of an M = 2 
U(N) gauge theory with L flavors and masses determined by the roots of B(z) (the lo- 
cations of the punctures). The curve is not generic: the further condition that f(z) be 
well-defined on E defined by y 2 = Q(z) implies that 

y 2 sw = P 2 {z) + B(z) = Q(z)H 2 (z) , (3.5) 

where H(z) is another polynomial. Plugging f(z) as given by ( |3.3j ) and ( |3.5|) in the 
expression (|3.1|) for the Strebel differential we get 

B(z)^pi(z)+B(z) 

= P ^ Q '^^T^ +2(-P'(z)H(z) +P(z)H'(z))y^dz. 
^Q{z)B(z) &{z) 



Comparing with (|2.2| ) we see that this expression has unwanted simple poles at the zeroes of 
Q(z), so these need to cancel for the differential to have the singularity structure of a Strebel 
differential. We impose this by requiring that Q(z) divides H(z), i.e. H(z) = Q(z)R(z) 
for some polynomial R(z), and obtain 



2ni ^dz = (-2P'(z)Q(z)R(z) + 3P(z)Q'(z)R(z) + 2P(z)Q(z)R ' (z)) dz . (3.7) 

Jd[ Z J 



Now note that, since ( p77|) was constructed as the logarithmic derivative of ( |3.3|) , its 
poles can only be simple poles at the distinct zeroes of B(z) and a simple pole at infinity 
with residue moo. Therefore we get the equation 



n-2 



[2P'(z)Q(z)R(z) - 3PQ'(z)R{z) - 2P(z)Q(z)R'(z)} = J\ (z - p^ 1 . (3.* 



i=0 



Substituting this into ( |3.7| ) we get the correct expression for the Strebel differential, 



dz = -^ dz. (3.9) 



If (L + rrioo) is not even, we multiply all residues by 2 and divide the differential by 2. 
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Note that the Strebel differential obtained this way has integral lengths. Let us see 
this in more detail. Define 6 to be 

VW) dz = ± dlog ( P p { *\ - y*" \ s l-4B. (3.10) 
2m \P(z)+yswJ 2tt 

Given z a and Zb to be any two zeroes of Q(z) (and therefore of ysw), we get 

~dz = ^-(e b -e a )ez, (3.ii) 

since both 9 a and 9b are equal to 2tvZ. So we have mapped the problem of finding Strebel 
differentials with integer lengths to that of solving the algebraic equation 



P 2 (z) + B(z) = Q 3 (z)R 2 (z). (3.12) 



3.2. Physical Interpretation 



Let us now comment on the physical interpretation of equation (|3.12|) . It is easiest to 
start with the curve ( |3?5| ) 



y 2 = P 2 {z) + B L (z) = Q 2n _ 4 (z)H 2 N _ n+2 (z) (3.13) 

where we have exhibited the degrees of the polynomials explicitly. This curve arises from 
a JV = 2 U (N) gauge theory with Nf massive flavors and tree level superpotential 

W tree = Tr W{§) + Qjm f f ($)Q f (3.14) 

where / and / run over the number of flavors Nf and 

n—2 _ i+i _ 

w & = E fti^ 1 ' m 'M = E m f,k zk ~ 1 ' ( 3 - 15 ) 

k=0 fc=l 

The degree of W(z) is n — 1 and is set by the fact that the supersymmetric vacua not 
lifted by ( |3.14j ) are those for which Q2n-4(z) = W'(z) 2 + f(z) where f(z) is a polynomial 
such that deg (/) = deg (W'(z))/2 — 1. These turn out to be N = 1 vacua. Finally, m(z) 
is a matrix of polynomials of size Nf x Nf. We have denoted the maximum degree of the 
polynomials in m(z) by I; lower degree ones are obtained by setting some of the fc 's to 
zero. 
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It turns out that the only information about the superpotential Qjm^{$)Qf which is 
relevant for the curve ( |3. 13 ) is the polynomial |3C,31,p^ 



B(z) = detmO). (3.16) 

Clearly, plenty of choices of m(z) can lead to the same B(z). Recall that we are interested 
in B(z)'b of the form Q 

n-2 

B(z) =az mo (z-l) mi Y[(z-Pj) mi . (3.17) 

i=2 

Two natural ways of obtaining such 5(z)'s are the following: 

• Nf = L = degi?(2;), i.e., Nf = 2N — m^ and m^(z) a constant diagonal mass matrix 
with mo masses equal to 0, mi masses equal to 1, and rrij masses equal to pj. 

• Nf = n — l and m^z) a diagonal matrix with polynomial entries z mo , (z — l) mi , and 
(z — pj) m j . 

The former leads to a theory with unbroken M = 2 supersymmetry if there is no 
W(z). Moreover, it has a large flavor symmetry classically. The latter, on the other hand, 
has a very small number of flavors and generically no special flavor symmetry. 

Now we are ready to consider ( |3.12| ). What we have done further is to tune the 
masses of the flavors (or in the second point of view the parameters of the polynomials) 
and the parameters of the superpotential to arrive at ( |3.12|) . Near any zero of Q(z) , the 



equation behaves as y 2 = x 3 , whose solutions are Argyres-Douglas points ]25j . These are 
the points where there are mutually non-local massless monopoles. Due to the presence of 
the superpotential the original monopoles condense and the new states give rise to M = 1 
superconformal theories in the IR. 

We stress that the factorization problem ( |3.12| ) is rigid, as it will be easily shown in 
the next section. From the point of view of the Strebel differential this means that, for a 
given set of integer residues, as we vary the location of the punctures, the differential has 
integer lengths only at isolated points in C n_3 . These points are the solutions of ( |3.12|) . 
In section 5 we will study the distribution of these points in the simplest nontrivial case 
of ?! = 4 and equal residues. 
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3.3. Alternative Construction 

Another equivalent way to arrive at ( |3.12| ) is to start with 

= ^dlog ( PW - Vf!W±^ (3.18) 

and realize that near a zero of order k of P 2 (z) + B(z) the abelian differential behaves as 
z k / 2 ~ 1 dz. It then follows that the corresponding quadratic differential behaves as z k ~ 2 dz 2 . 
Therefore, if we only want to have 2n — 4 simple zeroes and nothing else coming from the 
2N zeroes of P 2 {z) + B(z) , then k better be 3 for 2n — 4 of them and 2 for the rest. This 
immediately leads to (|3.12| ). 

The reason we took the longer route above was to make the relation to field theory 
manifest via ( |3.5|) and to lay down the basis for the construction of a differential equation, 
which will be the main tool for finding explicit solutions. 



4. Solving Argyres-Douglas Factorizations Using Differential Equations 

For the problem at hand, we have to find some polynomials P{z), Q(z), B(z) and 
R(z) satisfying the generalized Argyres-Douglas factorization problem ( |3.12| ). We choose 
P(z), Q(z) and R(z) to be monic. Let us count the number of unknowns and the number 
of equations. Imposing that P 2 (z) + B(z) factors as Q(z) 3 R(z) 2 is equivalent to requiring 
2n — 4 triple roots, which gives 2(2n — 4) conditions, and N — 3(n — 2) double roots, which 
gives N — 3(n — 2) conditions. The total number of conditions is thus N + n — 2. Now, the 
number of unknowns is given by N from P(z), n — 3 from p^-'s in B(z) and one more for 
a. This is exactly N + n — 2 , justifying our claim that the solutions are isolated points. 

Solving such polynomial equations is generically a tedious and complicated task. In 
this section we show that, by repeatedly solving linear equations, (|3.12|) can always be 
reduced to a small number (compared to 2N) of polynomial equations. We will look in 
particular at the case with four punctures, for which the number of final equations is two. 
Since the problem is rigid, i.e. there are only isolated solutions, the two equations are 
in two variables. A single equation for one of the two variables can then be obtained 
by computing the resultant of the two polynomials. This final equation factorizes into 
irreducible polynomials over Q , which have an interesting interpretation. 

Differentiation tricks are often useful in solving polynomial equations. A simple ex- 
ample is given by the Seiberg-Witten factorization corresponding to a maximally confining 
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point, where N — 1 mutually local monopoles become massless: P 2 (z) — 4 = (z 2 — 4)H 2 (z). 
(See appendix A.) 

Using the results of the previous section it is possible to find such a differential equation 
for P(z). Take two expressions for (p(z)dz 2 which are equivalent only when the factorization 
( P-12|) holds: for example, ( |3.6| ) and (|3.9| ). By comparing them, we get 



(P(z)B'(z)-2B(z)P'(z)) 2 mlQ(z 



(4.1) 



b{zY{p^z)+b{z)) n:=o 2 o^)' 

Substituting in this equation the definition of B(z) from ( |3.4j ) we finJ 

n — 2 r\TD( \ n ~ ^ 2 tl — 2 

[ P { z )Y^m i \[{z-p j )-2 d ^- Hiz-Pij) =ml Q(z)(P 2 (z)+aH(z-p j : 

i=0 j^i i=0 j=0 

(4.2) 

where we have set po = and pi = 1 as in previous sections. 

Recall from section 2.1 that imposing the residues at the poles pi to be equal to mi 
fixes Ti — l out of the 2n — 4 parameters of Q(z). Moreover, it turns out that the equations 
to be solved for the coefficients of P(z) are always linear and can be solved in terms of the 
unknown coefficients of Q(z) and B(z). We will discuss several examples in what follows. 

A particularly simple set of residues are those for which moo > mo + mi + . . . + m n _2 . 
In this case, P(z) can be found by solving the related differential equation 

n— 2 v n— 2 

/WE^H^-W)- 2 ^ H(z-p i ) = *Mzjm (4-3) 

and taking P(z) = [f(z)]+ (the polynomial part of f(z)). The reason is that when the 
condition moo > mo + wi + ■ • - + mn-2 is satisfied, the term proportional to a in (|4.2|) does 
not affect the highest + 1 powers of z. Therefore P(z) can be determined by dropping 
the a term completely and taking the square root. This will be used in Appendix B , 
where the series of residues (1,1,1, 3n) is solved for any value of n and for a particular 
location of the pole pi. 



7 This equation can also be derived starting from ( |3.8| ) if we use the factorization equation 
HD and the form of B(z) in 
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5. Examples 

In this section we consider some examples that illustrate how Strebel differentials with 
integer lengths can be constructed by solving the Argyres-Douglas factorization problem 
( p,12|) . Our main set of examples corresponds to the case of a sphere with four punctures 
but, as a warm up, we first look at the case with three punctures, which is trivial from the 
viewpoint of finding explicit Strebel differentials (as shown in section 2.2). 

5.1. Sphere With Three Punctures 

In section 2.2 we gave the explicit form of the Strebel differential on a sphere with 
three punctures and residues {77100,7710, mi}- In this section we want to show what the 
corresponding algebraic equation is. The factorization problem is 

P 2 (z)+B(z) = Q 3 2 (z)R 2 N _ 3 (z), (5.1) 

with 

B(z) =az mo (z-l) mi (5.2) 

and 2N = m^ + m + mi and L = deg B(z) = m + mi. 
Happily, the solution to the related algebraic equation 

P 2 {z) + G L (z) = (1 - z 2 ) H 2 N _ x {z) (5.3) 

was obtained in [^2j and is given in terms of linear combinations of Chebyshev polynomials. 
Here Gl(z) is a polynomial of degree equal to the number of effective flavors in the U(N) 
theory, as explained in Section 3.2. In the masses of all L flavors are arbitrary. 
The solution to ( |5"T5| ) is given by 

Nf 

Pn{z) = Y j v* cos((JV-i)0), 



i=0 
N f 

H N _i(z) = J2 



sm((N-i)6) 



sin (9 

i=0 

where z = cos 9 and the zVs are constants that will depend on the masses of the flavors 
or, equivalently, on the roots of Gl(z). 
We still need to impose that 

H N -i(z) = Q 2 (z) R N - 3 (z) 
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to get a solution of (O). This leads to two linear equations in the Vi variables, 



H N _ 1 (-1)=H N _ 1 (1) = 0, 

which can be solved easily. The only difficulty arises from requiring Gl(z) to have mo 
and mi coincident roots. The solution obtained in this way can be recast in the original 
form (|5.1| ) by an SX(2, C) transformation. This in principle leads to a full solution of the 
Strebel problem. 

Since the Strebel differential is already known for the case with three punctures and 
arbitrary residues and the factorization problem is already solved in the gauge theory, we 
move on to our first nontrivial set of examples. 

5.2. Sphere With Four Punctures 

There are very few explicit examples of non degenerate Strebel differentials known in 
the literature. To our knowledge there are only a few numerical results available for the 
case with four punctures, obtained setting all residues equal to one (this case is of interest 
for string field theory computations) p3|fj. In an explicit form for the differential is 
also given, for a very symmetric point in A-^0,4- 

Here, after a few comments about the most general case, we will study in detail the 
case of equal residues (n, n, n, n) , with bgN. For a given n we find Strebel differentials 
with integer lengths corresponding to all possible partitions of n in three natural numbers. 
When we rescale them to reinterpret them as differentials with residues (1, 1, 1, 1) we find 
that each solution (i.e. partition of n) corresponds to a different point of .Mo, 4- The 
picture that emerges from the examples that we work out explicitly is that, if all solutions 
were considered for all n, then the set of discrete solutions would fill up -Mo,4 , as expected 
from the discussion in section 2 . 

We will consider more examples in Appendix B. 

General Integral Residues: {moo, m ii m t} 

As explained in Section 3.1, the relevant factorization problem is 

P 2 N {z) + B(z) = Ql(z) R 2 N _ 6 (z) , (5.4) 



See [15 1 for very recent numerical results on the five point correlator in string field theory. 
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with iV = ~ (mo + mi + m t + m^) and 



B(z) =az mo (z-l) mi (z-t) mt . (5.5) 



Using the general form of the Strebel differential in ( |2.2|) and using the fact that the 
residues at 0,1, t are respectively mo, mi, mt, we can write (54(2) as 



Q 4 (z) = z 4 + Qiz 3 + Q2Z 2 + q^z + <? 4 , (5.6) 



with q 2 = 92(91, *), 93 = 93 (?!,*) and q 4 = 94(91, *) given by 



92(91,*) = 


9 

mf 










93(91, *) = 


/ mf 






94(91,*) = 


m o .2 

2 





(t - 1) + (^f + ^-(t-l))t-(l + q 1 +t + q 1 t + t 2 ) , 

+ (l + 9l) + ^(*-l)+fl-^V-4-( 1 +*)l *> ( 5 - 7 ) 



Therefore, the quadratic differential is specified by the two arbitrary complex parameters t 
and qi . Imposing that the differential is Strebel leads to two (real) independent conditions 
and thus determines q\ in terms of t. Since we are implicitly imposing the stronger condi- 
tion of integer lengths the complex parameter t is constrained to a discrete set of values, 
as we will see in the examples. 

Note that even after using the SL(2, C) symmetry to fix the location of three of the 
poles of <j)(z)dz 2 we still have a residual symmetry: we can exchange the positions of any 
two poles without changing the differential. As a consequence, 4>(z)dz 2 will be symmetric 
under the group of transformations generated by t — > 1 — t and t — > — 1/t . Moreover, t — > t* 
is also a symmetry of <p • These transformations have to be accompanied by corresponding 
transformations of q± and t so that the form ( |5.7|) remains the same. One can check that 
the following operations are symmetries: 

{z-^l-z, t^l-t, qi — >• -4 -qi}, 

{*--, t-7, 9i (5-8) 
z t t 

{z^z*, t->t*, qi^ql}- 

For every pair 

, 9i°' ) ) one can fi n d five more solutions to the factorization problem by 
performing the above operations. To mod out by these symmetries we will restrict t to the 
fundamental region bounded by the lines Re(£) = h, Im(t) = and the curve \t\ = 1. 
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Equal Residues: (n, n, n, n) 



We now restrict attention to the case mo = mi = = nit = n . As mentioned 
earlier, the simultaneous rescaling of all the residues does not affect the structure of the 
differential, therefore we can use this analysis to get information about the Strebel differ- 
ential with residues (1, 1, 1, 1) - after re-scaling by n - for different values of the complex 
structure parameter t . Although the Strebel lengths are integer valued for the (n, n, n, n) 
case, the rescaled lengths will typically not be integral. 

Let us first discuss some general features of the solutions we will find. Imposing that 
all the residues are equal leads to a simple form for the Strebel differential: 

9 • a A y^ 4 + qiz* - ( qi (l + t) - gg + qi tz + t 2 

2n i V a> dz = —n ; r- r dz . (5.9) 

vr z(z-l)(z-t) v ' 

A Strebel differential on a sphere with four punctures has a critical graph composed of six 
edges, with associated Strebel lengths. Once the four residues are specified, only two of 
the lengths are independent. A further simplification occurs when all residues are equal. 
The critical graph is a tetrahedron with all faces of equal perimeter and the lengths of the 
edges satisfy the relations (with notation that refers to Figure 1) 

i x = 4, £ 2 = £ 5 and £ 3 = £ 6 . (5.10) 

In what follows, we will label the Strebel differentials by three integers corresponding to 
the lengths 4,2,3, i- e - the lengths of the edges of a given face of the tetrahedron. 

In all the examples we considered we were able to find all solutions to the factorization 
problem 

P%{z) + a (z(z - l)(z - t)) n = Ql(z)R 2 N _ 6 (z) (5.11) 

with N = 2n and Q 4 (z) = z 4 + qtz s - (q^l + t) - 2t)z 2 + q t tz + t 2 . These are given by all 
possible sets of lengths consistent with residue n. In other words, all possible partitions 
of n into exactly three (strictly positive) integers. The number ps(n) of such partitions is 
given by the generating function 

3 00 

7 TT7 o TT = / PS (n)x n . (5-12) 

(x - l) 3 (x + l)(x 2 + X + 1) L^ LJ - iy - > V ) 

n=0 
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We obtained these solutions using the differential equation ( fO| ) , which for this case 
becomes 



(nP(z) (3z 2 - 2(1 + t)z + t)-2 z(z - l)(z - £)) 



n 2 Q 4 (z) (p 2 {z) + a (z(z - l)(z - t)) n ) 



(5.13) 



We expand both sides of the equation and compare the coefficients of the various powers of 
z . We find that the highest power z 2N+A always gives a trivial condition. Comparing the 
next iV coefficients we find iV linear equations for the coefficients of P{z) (this is always 
true because of the P'(z) term in the equation). The next condition, coming from the 
coefficient of z N+3 , gives a linear equation for a: this is not obvious, but a careful analysis 
shows this to be always true. At this point, all coefficients of P{z) and a are determined 
in terms of qi and t. 

Taking two more equations, say the coefficients of z N+2 and z N+1 , one discovers 
that they are nonlinear polynomial equations in two variables. The fact that they are 
nonlinear is actually good, since we know that we are supposed to find ps{n) solutions, 
with P3(n) given by ( |5.12j ). Taking the resultant of these two equations we find a single 
polynomial equation for t: let us call it Resi(t). The equation Resi(t) = gives us the 
solutions for t , but it also yields many spurious solutions. The reason is that there are 
more coefficients in ( |5.13| ) that will constrain the solution. Taking two more equations, say 



from the coefficients of z N and z N 1 , and computing a second resultant Res2(t) one can 
discard spurious solutions by computing 

fit) = GCD [Rest (t) , Res2 (£)] . 

The t values that solve the factorization problem are those that solve the equation f(t) = 0. 
In all our examples this was enough to discard spurious solutions, but in general one can 
continue this process by taking the coefficients of the lower exponents of z. 

We still need to compute q\. This is done by first computing the resultants Resi(q\) 
and Res2(qi) of fit) with the coefficient of z N+3 and z N+2 respectively (both of which 
depend on qi and t). The spurious solutions can be discarded by computing the GCD of 
these resultants: 

g( qi ) = GCB[Res 1 (q 1 ),Res 2 (q 1 )] . 
19 



As before, the q± values that solve the factorization problem are those that solve the 
equation g(qi) = 0. It turns out that both / and g are of the same degreei. 

We now turn to study particular cases. The cases n = 1, 2 are not interesting because 
from (|5.12| ) we see that P3(l) = and ^3(2) = 0. The way to see this from our factorization 
problem (|5.11] ) is that the degree of Pn[z) must be at least N = 6 . Therefore we need to 
take n > 3 . 

• n = 3 

For n = 3 the solution of the factorization problem 



is given by 



P§(z) +az 3 (z - l) s (z - t) 3 = Ql(z) 



1.^3 n . 2 , 64z /r ^ t , 

t = -+z-, qx = -2- l — and a = -^= (5.14) 



and the explicit form of the polynomials is 

P 3 (z) = 1 + (-1 + iV3)z + -(1 - iV3)z 2 + |(1 + iV3)z 4 + (-3 - iV3)z 5 + z 6 , 

9- 4-1 ( 5 - 15 ) 

Q 4 (s) = + (-2 - -^)z 3 + (1 + iV3)z 2 -2L z + -{-\ + iV3) . 

a/3 v 3 2 

The Strebel differential is 

0(^2 = _ 9 ^2 

47r2 z^{z-\y{z-\-i^Y 

The zeroes are located at the vertices of a regular tetrahedron with all the edges of length 
one in the Strebel metric. This is the most symmetric point in .Mo, 4 , the same found in 
||33j| , as mentioned at the beginning of the section0. 



9 This follows from the fact that each factorization problem has, associated to it, a unique 
number field over which the equation factorizes. This statement is related to the connection with 
Belyi maps and dessins d'enfants explained in section 7.2. 
10 Our conventions are such that —qi maps to the parameter a in [33]. 
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Figure 2: Location of poles and zeroes of the Strebel differential with residues (3, 3, 3, 3). 
• 4 < n < 9 

For 4 < n < 9 the values of (t,q±,a) and the corresponding lengths are listed in 
Table 1 . For all the values of n in the table the form of the polynomials P{z) and R(z) is 
easily found, but the expressions get cumbersome. Note that, whenever two of the three 
lengths are equal, the location of the third pole t is fixed to be in the line Re(i) = |. This 
automatically fixes the real part of qi to be Re((/i) = —2. 

Finally, let us make some remarks about something interesting that happens when n 
is not prime. Suppose that n = rs with r, s > 1 . Then from the uniqueness of the Strebel 
differential it must be that for each solution to the problems with n = r or n = s we get a 
solution to the problem with n = rs . This is indeed the case. Take for example n — 6 in 
Table 1 : the solution with lengths (2, 2, 2) is identical to the one for n = 3 with lengths 
(1, 1, 1). This observation seems surprising at first because the factorization problems for 
the two cases are quite different. However, this is a phenomenon already encountered in the 
study of Seiberg-Witten curves and is a consequence of the multiplication map introduced 
in H . 



Relation To Residues (1,1,1,1) 

As mentioned earlier, the solutions listed in Table 1 obtained by choosing integer 
residues (n, n, n, n) can be recast as solutions to the (1, 1, 1, 1) problem after a suitable re- 
scaling in each case. In Figure 3 we have plotted all values of t from Table 1, interpreting 
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the complex t plane as the moduli space .Mo, 4 with residues (1, 1, 1, 1) , a subspace of the 
full decorated moduli space -Mo,4 x • 

In Figure 3 , we have used the residual symmetry transformations (|5.8[) to bring those 
points into the fundamental region, whenever necessary (see section 5.2). It is clear that 
one can get an arbitrarily large number of points in this region by solving polynomial 
equations with higher n values. Note that each of these points is a well defined expansion 
point, around which a perturbative analysis similar to that performed in can be carried 
out. It would be interesting to use such an analysis to find, at least numerically, qi(t) as 
a function of the complex structure parameter t . 

We can see in this example an explicit realization of the general discussion of Section 
2.4 , where we showed that it is possible to approximate any Strebel differential, to arbi- 
trary accuracy, by solving the problem for integer lengths. Looking at Figure 3, one can 
observe that the points associated to lengths (l,r, s) follow a regular pattern, forming a 
grid parametrized by the integers (r, s). We have tried to make the grid more evident in 
Figure 4: only the points corresponding to the solutions that we have found explicitly are 
shown, but one should imagine an infinite grid. The same observation holds for the points 
associated to lengths (2,r, s) and we can presume it to hold more generally for lengths 
(p, r, s), for any integer p . 

Note that the lattice (2, r, s) contains the lattice (l,r,s) and it is finer. This has a 
simple explanation. Take the point corresponding to the lengths (1, 1, 2) : this is a node of 
the (1,7", s) lattice, but we know that the same value of t is guaranteed to give a solution 
for lengths (2,2,4) as well - or any other multiple. We conclude that the point (2,2,3) 
of the (2, r, s) lattice must lie between the origin and the first point of the (2, r, s) lattice. 
This simple argument shows how for increasing p the (p, r, s) lattice contains all lattices 
corresponding to p < p and it is finer. In the limit p — > oo we get a dense set of points 
covering the fundamental region of the moduli space. 
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n 


t 




a 


Lengths 


3 


- 4- 




64i 
3\/3 


(1,1,1) 


4 


2 + l 2v^ 


-2 -iV2 


-108 


(1,1,2) 


5 


1 + z 11 


2 z 2 f^ 


65536i 
2025\/T5 


(1,2,2) 




i + z 2.887948 


-2 - 1.527344 


1734.654912 z 


(1,1,3) 


6 


I _|_ jvf 
2^*2 




1024 
27 


(2,2,2) 




i + z 4.244172 


-2 - z2§ 


43595.583339 


(1,1,4) 




-0.127288 -z 0.622086 


-0.625270 + z 0.793701 


5.79167 + z 33.7770 


(1,2,3) 


7 


i +z 1.292013 


-2 - i - 1.314895 


-277.399584 % 


(2,2,3) 




0.220923 + z 0.294611 


-1.036986 -z 0.591846 


3.207583 + i 6.143871 


(1,2,4) 




i + z 0.341120 


-2 -z 0.697929 


-7.030850 z 


(1,3,3) 




| + z 5.841390 


-2 - z 1.623197 


-1.569485 x 10 6 


(1,1,5) 


8 




-2-i-\/2 


-2916 


(2,2,4) 




| + z 0.602725 


-2 -z 0.982567 


-23.372595 -z 3.803219 


(2,3,3) 




| + z 7.681492 


-2 - z 1.646274 


-7.639365 x 10 7 z 


(1,1,6) 




-0.106193 -z 0.309236 


-0.304343 + z 0.458200 


6.511956 + i 18.582505 


(1,2,5) 




-0.646518 - i 1.275342 


-0.617371 +i 0.981889 


-2202.314286 - % 810.756832 


(1,3,4) 


9 


2^*2 


- 2 -^73 


16384 i 
811\/3 


(3,3,3) 




i +z 2.298916 


-2 - z 1.480600 


41259.204415 z 


(2,2,5) 




i +z 0.271606 


-2 -z 0.592769 


0.000050 + z 6.353074 


(1,4,4) 




0.945257 + 0.754441 


-2.965934 -z 0.988093 


-107.573247 -z 19.721375 


(2,3,4) 




0.123306+ zO.189777 


-0.626671 -z 0.396566 


-4.130947 + z 4.589470 


(1,2,6) 




i +z 9.765341 


-2 - z 1.662029 


4.824601 x 10 9 


(1,1,7) 




1.969844 + i 1.659162 


-3.384879 -z 1.030218 


10737.225688 + i 29575.094989 


(1,3,5) 



Table 1: We list the t, q\ and a values for the (n,n,n,n) problem. The lengths £\, £2 and £3 of 
Figure 1 are specified in the last column: we get all partitions of n into exactly three integers. 
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n 
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n 


= 4 




n 


= 5 


• 


n 


= 6 


O 


n 


= 7 


® 


n 


= 8 


• 


n 


= 9 



(1,2,6) 



(1,2,5) 



(1,3,5) 



(1,1,7) 







(1,1,6) 



(1,1,5) 



(1,1,4) 



(1,1,3) 




Figure 3: Moduli space .Mo, 4 for residues (1, 1, 1, 1). The points in the plot are obtained solving 
the Strebel problem with residues (n, n, n, n) . The idea is to rescale by n to obtain a differential 
with the desired residues. Next to each point we have specified the corresponding Strebel lengths 
£i, £2, £3 (see the discussion around equation ( 5. JO| ) ). 



24 



6. Relation To Gauge Theory And String Theory 

In this section we expand upon the relation between Strebel differentials and M = 2 
gauge theories deformed by a tree level superpotential found in Section 3. We mention 
some of the problems on the gauge theory side that can be easily solved by using the 
relation to Strebel differentials. We also comment on some intriguing connections to string 
theory by geometrically engineering the gauge theories. 

6.1. The Generating Functional Of Chiral Operators T(z) 

Strebel differentials with integral lengths, as defined in ( |3.18| ) and (|3.3| ) , are intimately 



related to the abelian meromorphic differential T(z)dz defined in |p2[ . Using 

y 2 sw = P 2 (z)+B(z) 

we get 

2niyW)dz = dln( P ( z) - ysw ) 

v 1 ! \P(z) + yswJ 

= dlnB(z)-2dln(P(z)+y sw ) ^ 

= d\nB(z) -2T(z)dz. 

Since the first term does not contribute to a period integral unless it encloses a zero of 
B(z) (a pole of the Strebel differential), the factor of 2 in the relation between (p and T(z) 
ensures that the lengths of the Strebel differential are identical to the periods of T(z) . 



Li = — <p T(z) dz = Ni 
2ni J Ai 1 } 



for some choice of 1-cycles A{ that do not enclose any of the poles. As described in 
these integers Ni parametrize the classical vacua of a gauge theory in which the gauge 
group U(N) is broken from U(N) to n"=i U(Ni). 

It is important to note that it does not seem possible to interpret the full Strebel 



differential a/ <f>(z)dz as the generating functional of chiral operators for a physical theory. 
On the face of it, it would seem possible that our configuration with poles on the upper 
and lower sheet can be deformed (by moving the poles through the cuts) into a pseudo 
confining phase, discussed in |24|], in which all the poles are on the lower sheet. However, 
one can check that the starting point itself is unphysical, as some of the period integrals 
can be shown to be negative. In the gauge theory, these are interpreted as the rank of the 
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Figure 4: On top: Lattice of points with lengths (2, r, s) in Mo,a- The lattice (1, r, s) is contained 
in this and it is shown with solid lines. The dashed lines show the refinement of the lattice obtained 
including all the (2, r, s) points. The points marked by question marks are predictions. Lower 
left: Embedding in Z 2 of the lattice (1, r, s). Lower right: Embedding in Z 2 of the lattice (2, r, s). 
To compare the two lattices we need to shift the coordinates as (r — p,s — p) for the general case 
of the (p, r, s) lattice. We see clearly from the examples p = 1 and p = 2 that increasing p gives a 
rehnement of the lattice. 
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subgroups Ni into which the original U(N) is broken and so we do not attempt a direct 
interpretation of this sort. 

6.2. Counting Argyres- Douglas Points 

Consider the gauge theory described in section 3.2. This is a Af = 2 U(N) gauge 
theory with matter deformed by a tree-level superpotential of degree n — 1 for the adjoint 
scalar $ and a special superpotential Qjm^(^)Q^ for the coupling of flavors to <&. 

In section 3.2 we argued that the problem we would like to consider is the one that 
leads to a factorization of the form 

P 2 (z)+B(z) = Q 3 (z)H 2 (z) (6.2) 

with Q(z) = W'(z) 2 + f(z) where f(z) is a polynomial of degree n — 3 and W'(z) is the 
derivative of the superpotential function with degree n — 1. 

At these points in the Af = 2 moduli space something special happens: mutually 
non-local massless particles appear and therefore the Af = 2 theories are believed to be 
superconformal. After adding the superpotentials, supersymmetry is broken to Af = 1. 
The mutually local monopoles condense and the extra states might lead to Af = 1 super- 
conformal theories. The physics of this points might be complicated. However, we would 
like to concentrate on the problem of counting such points. 

The question is then, how many vacua are there with this property? 

Given just the problem of counting solutions to (|6.2| ) , finding the answer seems very 
hard. However, from the correspondence to Strebel differentials and metric ribbon graphs 
the problem can be easily converted into a combinatorial problem: to count all possible 
ribbon graphs with integer lengths determined by the given residues that can be drawn on 
a sphere! 

Let us consider some simple examples. Set n = 4. The factorization problem is 

P%(z) + az m °(z-l) mi (z-t) m * = Q 3 (z)R 2 (z). (6.3) 

The problem is to count the number of ribbon graphs with four faces, only trivalent 
vertices and such that the perimeter of the faces is given by {m^, mo, mi, to*} where 
nioo = 2N — m — ui\ — m t . The first step is to determine the possible topologies of 
the ribbon graphs. It turns out that in this case there are only five different topologies. 
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These are shown in Figure 5. One of them, the tetrahedron in Figure 5E , was already 
encountered in section 2.3. 

We now choose some particular families of {moo, mo, m i, m t} as examples. 

• {n, 7i, 7i, n}: This case is by now very familiar. The only possible topology is a 
tetrahedron and the number of ribbon graphs is given by the number of partitions of n 
into exactly three integers, i.e. 

oo 

^2p 3 (n)x n = x 3 + x 4 + 2x 5 + 3x 6 + 4x 7 + 5x 8 + 7x 9 + ... (6.4) 

n=0 

• {1, 1, 1, 2n + 3}: There is also only one possible topology in this case, that of Figure 
5A . It consists of three circles, each one connected by a line to a trivalent vertex. The 
circles must have circumference one. The number of possible ribbon graphs is given by the 
number of partitions of n in exactly three integers, i.e. p(n) . 

Examples for the other three possible topologies can also be easily constructed. 



6.3. Relation To String Theory 

The way we have written our Strebel differential is also closely related to the Seiberg- 
Witten differential of the corresponding M = 2 theory. For example, looking at equation 



12) of [35[] we see that the SW-differential can be written (after a suitable shift) as 



Xswdz = z din [ —J, 6.5 

\P(z) +ysw{z)J 

where y% w {z) = P 2 (z) + B(z) . 

From the work of we know that, geometrically engineering this theory in a type 
IIB superstring theory on a non-compact Calabi-Yau three-fold that undergoes a geometric 
transition, one finds that the SW-differential is given by 

Xswdz = z h(z)dz with h(z)dz = / H3 , (6.6) 

where H3 is the type IIB three-form field strength H^s + tHrr and r is the complex- 
ified string coupling. This means that for the geometries realizing the Argyres-Douglas 
singularities considered in Section 3, h(z)dz becomes a Strebel differential! 
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Figure 5: The five topologies of the ribbon graphs with four faces. 

This fact might have some connection with the counting of BPS states performed 
in |36]]. In that case finding BPS dyons in Argyres-Douglas superconformal theories was 
equivalent to finding special lagrangian cycles, whose definition contains the condition 
Im(e ia O) = , for some fixed real a and with O being the holomorphic three form of 
the Calabi-Yau three-fold. This is similar to the condition for a horizontal trajectory 
discussed in Section 2 . Thus the number of BPS states is roughly speaking the number of 
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horizontal trajectories. It is interesting to note that the condition for being Strebel could 
be interpreted as the condition to have a finite number of open horizontal trajectories and 
hence a finite number of BPS states. An important point to note is that the analysis of 
Shapere and Vafa is carried out by first zooming into the Argyres-Douglas point, which in 
our case would correspond to the SU(3) AD curve y 2 = x 3 , and then counting BPS states 
after deforming away from the AD point, by turning on relevant deformations so that the 
Seiberg-Witten curve becomes y 2 = x 3 + ux + v . By contrast, our differential is Strebel 
only at the AD-point and we retain information about the bigger U(N C ) gauge theory in 
which the AD point is embedded. Therefore it is not straightforward to identify objects 
between the two set-ups. It would be interesting to clarify the precise relationship between 
the two discussions. 



7. Strebel Differentials From Belyi Maps 



A surprising property of Strebel differentials was discussed in []27j : it was shown that 
a Strebel differential with integral lengths can always be constructed as the pullback by a 
rational function (3(z) of a meromorphic quadratic differential on a sphere with only three 
punctures: 

Here the map 

(3: E ^P 1 

must be a Belyi map [2(| , which means that it must satisfy the property of having exactly 
three critical values, at {0, 1, oo}. Let us see how this comes about in our construction. 
Consider the change of variables 

C= ^(l + cos#). (7.2) 



After this change of variables the differential on the r.h.s. of (|7.1| ) takes the form 



<P(z)dz 2 = ^de 2 (z). (7.3) 



If we equate this to our expression for the Strebel differential in ( |3.18| ) we get 



9 « W _ p ( z ) ~ VP 2 W+BU 



e% v [z) = _^ v w v L . (7.4) 

P{Z) + y/P 2 (z)+B( Z ) 
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Note that this is exactly the change of variables performed in ( 3.10 ). Substituting this in 
the expression for (|7.2|) , we get an explicit expression for our candidate Belyi map in terms 
of the polynomials that solve the factorization problem (|5.4j ) : 

Q 3 (z)R 2 (z) P 2 (z) 
m ~ B(z) ~ 1 + B(z) ■ {7 - b) 

Let us now show that the map ( |7.5| ) satisfies the right properties: 

• From the two equalities, it is clear that the critical points of (3 are the zeroes of 
Q(z), R(z) and P(z). 

• The critical value at the zero of Q(z) or R(z) is zero, while the critical value at a 
zero of P(z) is equal to 1. 

• oo is a ramification point on P 1 and the pre-image of oo is given by zeroes of B(z). 
The map ( |7.5|) is ramified on P 1 at exactly the points {0, 1, oo} and is thus a Belyi map. 

7.1. Example : (2,2,2) 

As an example, let us apply this formula to the simple case of the sphere with three 
punctures, choosing residues (2,2,2). This case is discussed in |2~T|] . The polynomials now 
satisfy the equation 

Pi(z) + az 2 (z-l) 2 = QUz). 
Solving the differential equation (|4.2| ) for this case, we find the unique solution 

P 3 (z) = 1 - ^z - ^z 2 + z 3 Q 2 (z) = l-z + z 2 and 
and substituting this into (|7.5| ) we get 

4 (z 2 -z + l) 3 



(3(z) 



27 z 2 (z-l) 



2 ' 



which is the answer quoted in E7 



7.2. Some Comments on Belyi Maps, Children's Drawings and Seiberg-Witten Theory 

It turns out that the surprising connection we have found between Belyi maps and 
Argyres-Douglas curves (in ( ^.12| )) can be generalized to any Seiberg-Witten curve that 
develops an isolated singularity. The precise connection will be explored in a forthcoming 
publication [[28 . Here we just collect a few well known facts about Belyi maps and comment 
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briefly on the connection. The proofs of most of the technical statements to follow can be 
found in IR7 . 



If we are given a Belyi map in terms of rational functions, say 

= (7 ' 6) 

it is always possible to associate to it a diagram on the Riemann sphere. As explained in 
|Tj],^7| , these graphs (usually referred to as Grothendieck's Dessins d'Enfants or "children's 
drawings") are given by the inverse image under (3 of the interval [0, 1] on P 1 : 

D = p- 1 ([0,l]). (7.7) 

A dessin is said to be "clean" if all the ramification indices at the pre-images of 1 are 
exactly equal to two. We see that this is satisfied by ([7.5|) and so the dessins we construct 
using ( |7.7p are clean. More generally, the Belyi map in (|7.6| ) leads to clean dessins iff the 
polynomials in (|7.6| ) satisfy the equation 

A(z)-B(z) = P 2 (z), (7.8) 

for some P{z). It can be shown that a clean dessin has as many vertices as there are distinct 
zeroes of the Belyi map and as many open cells as there are distinct poles. The valency of 
each vertex is given by the order of the corresponding zero. Similarly, the valency of each 
open cell (i.e. the number of edges that bound it) is given by the order of the pole. 

It is not difficult to see from these definitions that, when the polynomial equation ( |7.8| ) 
describes the Argyres-Douglas singularity in ( |3.12|) , the dessin defined by (|7.7| ) coincides 
with the critical graph of the Strebel differential. More generally, we can identify any 
polynomial equation ( |7.8| ) that gives rise to a clean Belyi map with a Seiberg-Witten 



curve. From (|7.7|) above, there will be a children's drawing associated with any such 
Seiberg-Witten curve. An example that will be explored in detail in this direction [[28] is 
that of pure U(N) gauge theory, for which B(z) = — 4A 2Ar , with A the strong coupling 
scale of the theory. 
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7.3. Assigning Lengths To Drawings 

From our discussion at the end of Section 3.1, it is clear that the Strebel differentials 
coming from Belyi maps have integral lengths. We now extend this by associating a general 
differential ^ to any drawing D that can be obtained from a Belyi map using (|7.7|) . Using 
4>d one can assign lengths to the edges of the drawing. Let us see how this comes about. 

It follows from the definition of the Belyi map in (|7]6|) that every edge which goes 
between any two successive zeroes of A(z) (pre-image of 0) has on it a pre-image of 1, 
which corresponds to a zero of P(z). This is shown below in Figure 6 for the specific case 
A(z) = Q\(z), B(z) = B^(z) and P(z) a degree 6 polynomial: a case discussed in Section 
5 as the (3, 3, 3, 3) example. 




O zeroes of P(z) 
# zeroes of A(z) 



Figure 6: Location of the pre-images of and 1 under the Belyi map . Pre-images of correspond 
to zeroes of A(z) and pre-images of 1 correspond to zeroes of P(z). 

Let us use the definitions in (|7.3j ) and ( |7.4| ) for polynomials P(z) and B(z) that satisfy 
([7.8|) to construct the differential 4>d '■ 

From O and O it follows that 



>b 



^- [ y/<P D (z)dz= -!- / d0, (7.10) 



2ni J a v r v ' 2tt 

where a and b are the zeroes of A(z) and P(z) respectively and 9 a and 9^ are the corre- 
sponding ^-coordinates. From ( [7.9|) we see that at a zero of A(z), we get 9 a = mod 2n 

33 



and at a zero of P(z), we get 0& = tt mod 2%. Substituting these values into ( |7.10| ) leads 
to 

-i- / V^5)dz=i + Z. (7.11) 
As shown in Figure 6, when a and 6 are adjacent on the graph, we get 

8. Concluding Remarks 

We have shown that finding Strebel differentials with integer lengths is equivalent 
to solving an algebraic problem. This algebraic problem is the same as that of finding 
generalized Argyres-Douglas singularities in the Coulomb moduli space of an Af = 2 gauge 
theory, which corresponds to a factorization of the Seiberg-Witten curve of the form ( |3.12|) . 

This correspondence between Strebel differentials and M = 2 gauge theory turns 
out to be quite useful: the relation to Strebel differentials allowed us to write down the 
differential equation (|4.2|) that was instrumental in solving the factorization problem. 

We showed, using the correspondence of Strebel differentials to metric ribbon graphs, 
that any Strebel differential can be approximated by one with integer lengths divided by 
an appropriate integer. In particular, for the example with four punctures, by solving the 
(n, n, n, n) case with integer lengths we could obtain the set of points in the moduli space 
with redidues (1, 1, 1, 1) that lead to rational Strebel lengths. One can see from Figure 4 
that the lattice of points that solve the (n, n, n, n) problem get finer and finer as one solves 
for higher and higher n values. Taking n to be very large, one can intuitively understand 
how the solutions the Strebel problem with rational lengths form a dense subset of the 
moduli space. It would be interesting to develop a method for finding an interpolating 
function q± = qi (t) from the data we have. 

The correspondence between Strebel differentials and metric ribbon graphs also al- 
lowed us to find the number of solutions to the factorization problem ( |3. 12|) , something 
that is not at all obvious from the gauge theory side. For all residues equal to n, we found 
a particularly simple answer: there are as many solutions as there are partitions of n into 
exactly three strictly positive integers. 

The Strebel differential also has a simple interpretation in the M = 2 gauge theory: 
it is related to the generating function of scalar correlation functions T(z) dz as shown 
in ( |6.1|) . It has an even more striking interpretation in string theory, in the context of 
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geometric engineering: it is the if-field integrated over an S 2 in the non-compact Calabi- 
Yau. It would be extremely interesting to explore this in more detail. In particular, we 
found that the conditions imposed on a general quadratic differential to be Strebel were 
very similar to the conditions found in |3(| to find special lagrangian three manifolds in 



a non-compact manifold. Interestingly, the Calabi-Yau in question is conjectured to be 
holographically dual to Argyres-Douglas superconformal theories 

Finally, from the discussion in Section 7, it is clear that the relationship between 
Seiberg-Witten theories and clean Belyi maps is more general and not restricted to the 
special class of Argyres-Douglas singularities discussed in this paper. The relation between 
Belyi maps and the children's drawings raises the interesting question of understanding the 
role of such drawings in gauge theory. In this note we have been content with describing 
how our understanding of the relationship between Argyres-Douglas curves and Strebel 
differentials leads to a differential equation technique to solve the problem of factorizing 
the Seiberg-Witten curve as in ( ^.12| ). We will explore the more intriguing relationship 
between children's drawings and Seiberg-Witten theory in a forthcoming publication ||28|| . 
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Appendix A. Maximally Confining Vacua 

The maximally confining vacua P,|4C)[1 are isolated singular points of the Seiberg-Witten 



curve of pure Af = 2 SU (N) SYM which correspond to the factorization problem 

P 2 {z)-A={z 2 -A)H%_ x {z). (A.l) 

In this appendix we show how differentiation gives a simple way of finding the solutions. 
After differentiating once on both sides we get 

2P N (z)P' N {z) = H N - X (z) [2zH N _ x (z) + 2(z 2 - 4)H' N _ 1 (z)] . (A.2) 
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Note that a root of Pn( z ) on the l.h.s. cannot possibly be a root of Hn-i(z) on the r.h.s., 
because that would contradict ( |A.1|) . Therefore, all roots of Hn-x(z) must be roots of 
P' N (z). Since the two polynomials Hjq-i{z) and Pn{z) have the same degree and are both 
monic, it must be that 

P' N (z) = NH N - 1 {z). (A.3) 



Taking this into account, equation (|A.2|) implies 



NP N (z) = zH N - X {z) + (z 2 - A)H'n-i(z) (A.4) 
and using (|A.3|) to solve for Hjq-\{z) we get 

N 2 P N (z) = z^(z) + (, 2 - 4)^g=i(z) • (A.5) 
This differential equation is easily recognized as the Chebyshev equation and it has two 



independent solutions T^(z) and \/ z 2 — AUn-i{z) , where Tjy(z) and Un-x(z) are the 
Chebyshev polynomials of first and second kind respectively. Since Pn(z) is a monic 
polynomial we conclude that 

P N {z) = 2T N (|) and H N _ x {z) = U N _ X (|) . 



Appendix B. Strebel Differential On A Sphere With Residues (l,l,l,3n) 

For this case one can solve the factorization problem ( |5.4| ) for any odd value of n . 
(For even n, one cannot define a polynomial factorization problem as the degree of N is 
no longer integer.) Since we require the lengths to be integral, the only possible way to 
draw the critical graph -such that the sum of positive integers is equal to 1 for three of the 
loops - is to have three of the edges be just curves that encircle the location of the poles 
at (0, l,t) . Since all vertices are trivalent, this fixes the form of the critical graph to be 
Figure 7 . The residue at infinity being 3n then leads to the condition 

5> = ^-i). 

As we argued in the main text, the number of solutions (distinct values of t and q±) of 
this problem is P3(f [n — 1)), where P3(n) is the partition of n into three non-zero positive 
integers. However, it is very easy to find one solution of the factorization problem: if we 



36 



require that all Zj be equal, then the points (0, 1, t) form an equilateral triangle, which fixes 
t - and hence qi - to be 



1 y/3 
- + i 

2 2 



Qi 
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(B.l) 



This completely determines in ( |5.7| ). Moreover, since this case satisfies the condition 
that > m>o +mi +m t , one can use ( |4.3| ) and solve for Pn(z) completely. Substituting 
Pn(z) and Q±(z) into the full differential equation (|4.2| ) one can solve for a. Plugging all 
this back into the original factorization problem, we can find R(z). 




Figure 7: Location of poles and zeroes of the Strebel differential with residues (15, 1, 1, 1) . The 
dashed line around the critical graph represents the contour that computes the residue at infinity. 

In practice, this program is not easy to carry out. However, by solving the equations 
for the first few cases, we find that the factorization problem 



Pl {ri+1) (z) +az(z- l)(z -t)= Ql(z) R 



f(n-3)» 



is solved by the ansatz (p.l|) if a is given by the closed expression 



/ n+1 \n+l( n—1 \n—l 



-1) — 4i 



3 3 2™ n 2n 
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The Strebel differential for this problem has the form 

47T2 Z 2 {z _ mz _ 

where 

Q 4 (z) =z* + (-2 - ^=)z* + (1 + i V3)z 2 - ^- 2 {l + 3n*)z + JL(-1 + zv^) . 
The zeroes of Q±{z) are located at 

Zl = 1(3 + zv^) , z 2 = 1(3 + i>/3) + (3i + V3){n 4 - n 6 )? 



zs = 1(3 + iV3) - -^=(n 4 - n 6 )^, z 4 = 1(3 + zv^) + ^J^ ) 



One can explicitly compute the Strebel lengths and check that they are for |z i|, 
i = 1,2, 3. The notation is as in Figure 7 , where we have shown the critical graph for the 
particular case n = 5 . For this the edges have lengths zoi = 2 and za = 1 so the residue 
at infinity is (correctly) given by m^o = 15. 
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